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A new measurement of the branching ratio, -R e / M = r(-7r + —> e + u + n + e + v''y)/T('K + —► 
p + u + 7r + —> p + u 7), resulted in = (1.2344 ± 0.0023(sfat) ± 0.0019(sysf)) x 10 -4 . This is in 
agreement with the standard model prediction and improves the test of electron-muon universality 
to the level of 0.1 %. 

PACS numbers: 13.20.Cz, 14.40.Be, 14.60.St, 14.80.-j 


The standard model (SM) assumes equal electro-weak 
couplings of the three lepton generations, a hypothesis 
known as lepton universality which is studied in high 
precision measurements of 7and W decays. A 
recent measurement of B + —» K + l + l~ decays [l|, where 
l represents e or /x, hinted at a possible violation of e-/x 
universality in second order weak interactions that in¬ 
volve neutral and charged currents. The branching ratio 
of pion decays, i? e / M = r(7r —>• ez/(7))/r(7r —» fxv( 7)), 
where (7) indicates inclusion of associated radiative de¬ 
cays, has been calculated in the SM with extraordi¬ 
nary precision to be = (1.2352 ± 0.0002) x 10 -4 

Hi- Comparison with the latest experimental values, 
R e e x ^ = (1.2265 ± 0.0034(staf) ± 0.0044(sysf)) x 10" 4 0 
and R ex P = (1.2346±0.0035(siai)±0.0036(syst)) x 10“ 4 
[I], has provided one of the best tests of e-fx universality 
in weak interactions for the charged current, at the 0.2 % 
level giving sensitivity to new physics beyond the SM up 
to mass scales of 0(500) TeV[3(. Examples of new physics 
probed include R-parity violating SUSY [dj, extra leptons 
Q and leptoquarks Q. In this paper, we present the first 
results from the PIENU experiment, which improve on 
the precision of A)) 1 ) and the test of e-fx universality. 

The branching ratio R e / M is obtained from the ratio 
of positron yields from the 7r + —» e + v , ( , y) decay (total 
positron energy E e + = 69.8 MeV) and the 7r + —> fx + v( 7) 
decay followed by the fx + —> e + uV('j) decay (7T + —> p. + —> 
e + , E e + = 0.5 — 52.8 MeV) using pions at rest. FigureQ] 


shows a schematic view of the apparatus [§j in which a 
75-MeV/c 7r + beam from the TRIUMF M13 channel [l(| 
was degraded by two thin plastic scintillators B1 and B2 
and stopped in an 8-mm thick scintillator target (B3) at 
a rate of 5 x 10 4 7r + /s. Pion tracking was provided by 
wire chambers (WC1 and WC2) at the exit of the beam 
line and two ( x,y ) sets of single-sided 0.3-mm thick planes 
of silicon strip detectors, SI and S2, located immediately 
upstream of B3. 
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FIG. 1: Top half cross-section of the PIENU detector. The 
cylindrical Nal(Tf') crystal is surrounded by a cylindrical ar¬ 
ray of Csl crystals as described in the text. 


The positron calorimeter, 19 radiation lengths (r.l.) 
thick, placed on the beam axis consisted of a 48-cm (dia.) 
x 48-cm (length) single-crystal Nal(Tl') detector |ll| pre- 
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ceded by two thin plastic scintillators (T1 and T2). Two 
concentric layers of pure Csl crystals [12] (9 r.l. radi¬ 
ally, 97 crystals total) surrounded the Nal(Tt') crystal to 
capture electromagnetic showers. Positron tracking was 
done by an ( x , y ) pair of Si-strip detectors (S3) and wire 
chambers (WC3) in front of the NaI(T£) crystal. 

A positron signal, defined by a T1 and T2 coincidence, 
occurring in a time window 300 to 540 ns with respect 
to the incoming pion was the basis of the main trigger 
logic. This was prescaled by a factor of 16 to form an 
unbiased trigger (Prescaled-trigger). Events in an early 
time window 6 to 46 ns and events with E e + > 46 MeV in 
the calorimeter provided other triggers (Early- and HE- 
triggers), which included most 7r + —> e + v decays. The 
typical trigger rate (including monitor triggers) was 600 
Hz. 
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FIG. 2: Energy spectra of positrons in the time region 5 
to 35 ns without and with (shaded) background-suppression 
cuts (see the text). The vertical line at 52 MeV indicates the 
E cu t position. 

Events originating from stopped pions were selected 
based on their energy losses in B1 and B2. Any events 
with extra activity in the beam and positron counters 
(Bl, B2, T1 and T2) in the time region of 7 to 1.5 /xs 
with respect to the pion stop were rejected. About 40 % 
of events survived the cuts. A fiducial cut for positrons 
entering the NaI(T£) detector required a track at WC3 
to be within 60 mm of the beam axis to reduce electro¬ 
magnetic shower leakage from the crystal. 

The summed NaI(T£) and Csl energy for positrons in 
the time region 5 to 35 ns is shown in Fig. [2j The time 
spectra for events in the low- and high-energy regions 
separated at E cut = 52 MeV are shown in Fig. [3] Events 
satisfying the Early-trigger or Prescaled-trigger filled the 
low-energy histogram (Fig. [3ji) and HE-trigger events 
filled the high-energy histogram (Fig. [3]o) . There were 
4 x 10 5 7 r + —> e + u events at this stage. The raw branch¬ 
ing ratio was determined from the simultaneous fit of 
these timing distributions. To reduce possible bias, the 
raw branching ratio was shifted (“blinded”) by a hidden 


random value within 1 %. Prior to unblinding, all cuts 
and corrections were determined and the stability of the 
result against variations of each cut was reflected in the 
systematic uncertainty estimate. 

In the low-energy time spectrum, the main components 
were 7r + —> y + —> e + decays at rest (LI), /i + —> e + vv 
decays (L2, about 1 % of LI) after decays-in-flight of pi¬ 
ons (7 tDIF), and decays coming from previously stopped 
(“old”) muons remaining in the target area (L3): 

LI: F hl = e“ A V - e" A -‘) for t > 0, 

L2: Fl 2 = A M e -A V for t > 0, and 

L3: Fl 3 = A M e -A V for any t. 

The distribution coming from the presence of plural 
muons in the target area was estimated to be <0.01 %, 
and was ignored in the fit. The low-energy fraction of 
7r + —> e + v events due to shower leakage and radiative 
decays was also negligible in the low-energy time spec¬ 
trum fit. 
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FIG. 3: (Color online) Time spectra of positrons (thin line 
histograms) in the (a) low- and (b) high-energy regions sep¬ 
arated at E cu t. The notches at t = 0 ns are due to a veto 
for prompt pion decays, and the peak at -3 ns in (b) is due 
to positrons in the beam. Each curve, labeled with the corre¬ 
sponding component described in the text, indicates the am¬ 
plitude in the fit. LI and part of L3 significantly overlap with 
the data. The thick solid line in (b) for t < 0 ns shows the fit. 
The fit for the other regions is almost indistinguishable from 
the data and is omitted here. 

The primary time distribution component in the high- 
energy region was the 7r+ —> e + v decay (HI: Fhi = 
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A 7 r e _A,rt for t > 0). The amplitude of HI also included 
the high-energy portion (E e + > E cut ) of decay-in-flight of 
muons (/rDIF) following 7r + —> ul + v decay at rest, which 
was estimated by simulation |l3( to be (2.07 ± 0.06) x 
1CT 7 of LI. The major backgrounds (H2) in the high- 
energy region came from muon decays due to the energy 
resolution of the detector, radiative muon decays in which 
the 7 -ray raised the observed calorimeter energy above 
E C ut, and extra hits (pile-up) in the calorimeter with a 
flat time distribution ( e.g. due to neutrons from the pion 
production target). The H2 component had an identical 
time dependence to the low-energy spectrum (L1+L2). 
The contribution from L3 via the same mechanism was 
separately treated as a muon decay component (H3) to 
include other contributions of “old” muons. 

Radiative pion decays 7 r + —> fi + v 7 (branching frac¬ 
tion, 2 x 10 -4 0 ) followed by fi + —> e + vv decays could 
contribute to the high-energy region if the 7 -ray hit the 
calorimeter. The contribution of the extra 7 -ray to the 
observed positron energy varied with the time difference 
of the two decays. This contribution (H4) was simulated 
using the observed pulse shapes of the Nal(Tf) and Csl 
signals, and is shown by the dashed line in Fig. [HJo- The 
amplitude of this component was (4.9 ±1.0) x 10 ~ 7 of LI 
in the fit. 

The background in the region t < 0 ns was due to 
events with time distribution H3 and those in which a 
positron from an “old” muon hit T 1 in coincidence with 
a positron, from 7r + —> /r + —>• e + decay of the stopped 
pion, that missed T1 but hit the calorimeter, raising the 
observed energy above E cut . The shape of this time 
spectrum (H5), including the inverse combination, was 
generated by simulation using the observed pulse shapes 
and the energy distributions for the corresponding event 
topologies. The shape and the relative amplitude of H5 
are shown by the dotted line in Fig. [3 Jd. 

A pile-up cut based on the T1 waveform rejected events 
with two hits. However, events with two T1 hits within 
the double pulse resolution of T1 (AT = 15.7 ± 0.3 
ns) were accepted, and the probability for the measured 
positron energy to be E e + > E cut was high in those 
events. By artificially increasing the double pulse res¬ 
olution up to 200 ns, the amplitude of this component 
(H 6 ) was obtained and fixed to the “old” muon back¬ 
ground L3 in the fit. The H 6 component is shown by the 
full line in Fig. [3jy the uncertainty in R e / fl was 0.01 %. 

The free parameters in the fit for the low-energy region 
were the amplitudes of LI, L2, and L3. The time origin 
1 0 , which was determined using prompt events, was fixed 
in the fit. The choice of to did not affect the branching 
ratio as long as the amplitude of L2 was a free parame¬ 
ter. The free parameters for the high-energy region were 
the amplitudes of HI, H2, H3 and H5. The total y 2 of 
the high-energy and low-energy fit was minimized with 
a common to- The fitting region was -290 ns to 520 ns 
excluding the prompt region of -19 to 4 ns. 


The overall fit result is almost indistinguishable from 
the data and is not displayed in Fig. [3j except in (b) for 
t < 0 ns (thick solid line). No structure was evident in 
the plot of residuals of the fit. The raw branching ra¬ 
tio and its statistical and systematic uncertainties were 
= (1.1972 ± 0.0022(staf) ± 0.0005(syst)) x 1CT 4 
with y 2 /DOF = 1.02 (DOF=673). The systematic un¬ 
certainty includes uncertainties of the parameters and 
shapes in the fit and of small components excluded from 
the standard fitting function as listed in Table HI The 
branching ratio was stable for the fits with free pion and 
muon lifetimes, which were consistent with the current 
values 0. 

Some corrections applied to the raw branching ratio 
relied on simulation [131 ]. Pions were generated 0.5 m 
upstream of the detector according to the measured pion 
beam distribution. Small energy-dependent effects in the 
energy-loss processes of positrons change the relative ac¬ 
ceptances of low- and high-energy events. The ratio of 
the acceptances of the 7r+ —> [i + —> e + and 7r + —»• e + v 
decays was found to be 0.9991 ± 0.0003(s?/sf) for a WC3 
radius cut r <60 mm. 


The largest correction to the raw branching ratio was 
for the 7 r + —>• e + u events below E cut , which primar¬ 
ily arose from the response function of the calorime¬ 
ter. Because of structure in the response function due 
to hadronic interactions [16], which was not well repro¬ 
duced by the simulation, empirical measurements were 
performed. Special data, using a simplified setup consist¬ 
ing of T2 and WC1-3, taken with a 70-MeV/c positron 
beam at various entrance angles were used to deter¬ 
mine the response function. In order to obtain the 
fraction of the 7r + —>• e + v events below E cut for the 
full setup, the difference in the detector geometry, the 
7 r + —>• c + i/ angular distributions, and radiative pion 
decays were estimated using simulation. The fraction 
of the events below E cut = 52 MeV was found to be 
3.19 ± 0.03 (stat) ± 0.08(s?/st) %. Since a small contribu¬ 
tion to the observed fraction from low-energy positrons 
in the beam could not be ruled out, the tail correction 
obtained in this way was treated as an upper bound. 

In order to estimate the lower bound to the tail frac¬ 
tion, 7r + —>• /i + —>- e + events were suppressed using 
an early decay-time region 5-35 ns, pulse shape and 
total energy in B3, and measurements of the straight¬ 
ness of the pion track [175 ■ The resulting background- 
suppressed positron energy spectrum is shown by the 
shaded histogram in Fig. [2] The remaining background 
was subtracted from the spectrum using the fact that the 
background-suppressed spectrum in a low-energy region 
contained a negligible 7r + —>• e + v tail contribution. The 
area of the low-energy region was scaled to the full region 
(< E C ut) using the known background distribution. This 
resulted in a lower bound of 1.48±0.07(stat)±0.08(sj/st) 
%. Since the total energy cut used in the suppression 
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method tended to remove 7 r + —> e + v events with Bhabha 
scattering which resulted in larger energy deposit in B3, 
a correction of 1.48±0.02(s?/st) % obtained by simulation 
was added to the tail correction. Thus, the lower bound 
was 2.95 ± 0.07(stat) ± 0.08 (syst) %. Combining the up¬ 
per and lower bounds, a multiplicative tail correction of 
1.0316 ± 0.0012 was obtained. 

Possible energy-dependent effects on t Q were studied 
using positrons in the beam at momenta 10-70 MeV/c, 
and with positrons from muons stopped at the center of 
B3 by lowering the beam momentum to 62 MeV/c. The 
multiplicative correction from this effect was 1.0004 ± 
0.0005. Other uncertainties included are for possible trig¬ 
ger inefficiencies (±0.0003) and distortions due to pile-up 
and other cuts (±0.0005). 

Stability of the measured branching ratio was further 
tested for dependence on many parameters, such as fit¬ 
ting ranges, fiducial cuts, pile-up cuts and E cut , which 
provided confidence in the validity of the background 
functions and corrections. Figure |4] shows the depen¬ 
dence on E cut . The drop below 50.5 MeV is primarily 
due to the energy threshold of the HE-trigger. 



FIG. 4: Dependence of the branching ratio on E cu t with re¬ 
spect to the value at 52 MeV. The error bars indicate addi¬ 
tional statistical and systematic uncertainties. The variations 
indicated here are small compared to the statistical uncer¬ 
tainty of 23 x 10 -8 at E cll t = 52 MeV. 

Table Q] shows a summary of the fit uncertainties and 
corrections after “unblinding”. The measured branching 
ratio is tffff = (1.2344 ± 0.0023(s±rf) ± 0.0019(syst)) x 
10 -4 , consistent with previous work and the SM predic¬ 
tion. The present result improves the test of e-p, univer¬ 
sality compared to previous experiments by a factor of 
two: g e /gfi = 0.9996 ± 0.0012 for the charged current. 
Results using an order of magnitude more data and pos¬ 
sibly improved systematic uncertainty estimates will be 
forthcoming. 

This measurement also results in improved 90 % 
confidence-level limits 0 on the neutrino mixing param¬ 
eter U e i between the weak electron-neutrino eigenstate 


and a hypothetical mass eigenstate rn Vi [17f, \U e i \ 2 < 


0.0033/(p e — 1) in the mass region < 55 MeV, where p, 
is a kinetic factor found in Refs. 00. 



Values 

Uncertainties 



Stat 

Syst 

Rf/T ( 10 - 4 ) 

1.1972 

0.0022 

0.0005 

7 r,/i lifetimes 



0.0001 

Other parameters 



0.0003 

Excluded components 



0.0005 

Corrections 




Acceptance 

0.9991 


0.0003 

Low-energy tail 

1.0316 


0.0012 

Other 

1.0004 


0.0008 

Rf/T ( 10 - 4 ) 

1.2344 

0.0023 

0.0019 


TABLE I: The table includes the raw branching ratio with 
its statistical and systematic uncertainties, the multiplicative 
corrections with their errors, and the result after applying 
corrections. 
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